Hawking radiation of nonsingular black holes in two dimensions 
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In this letter we study the process of Hawking radiation of a black hole assuming the existence of a 
limiting physical curvature scale. The particular model is constructed using the Limiting Curvature 
Hypothesis (LCH) and in the context of two-dimensional dilaton gravity. The black hole solution 
exhibits properties of the standard Schwarzschild solution at large values of the radial coordinate. 
However, near the center, the black hole is nonsingular and the metric becomes that of de Sitter 
spacetime. The Hawking temperature is calculated using the method of complex paths. We find that 
such black holes radiate eternally and never completely evaporate. The final state is an eternally 
radiating relic, near the fundamental scale, which should make a viable dark matter candidate. We 
briefly comment on the black hole information loss problem and the production of such black holes 
in collider experiments. 



PACS number(s): 04.60.-m; 04.70.Dy; 98.80.Cq. 
I. INTRODUCTION 

Two of the outstanding questions in physics today are 
how should Einstein's theory of gravity be modified in 
regions of spacetime with high curvatures, and what is 
the final state of an evaporating black hole. The physics 
required to answer both of these questions rests in an as 
of yet unknown theory of quantum gravity. 

In general, theories which report to describe features 
of quantum gravity (including some string theories) con- 
tain higher derivative corrections to the Einstein- Hilbert 
action [1] - [6]. For example, perturbative string theory 
predicts an infinite series of correction terms to the Ein- 
stein equations 
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= 0. (1.1) 



These correction terms become important at energies 
near the string scale yV = M~ l . It is possible that such 
terms will drastically affect the structure of singularities 
in both cosmological and black hole spacetimes. Perhaps 
quantum gravity will some how remove the singularities, 
leading to geodesically complete spacetimes. 

Two examples of higher-derivative gravitational theo- 
ries which eliminate the singularity in the Pre-Big-Bang 
model (allowing for a successful graceful exit) are con- 
structed in [7,8]. These models utilize the Limiting Cur- 
vature Hypothesis (LCH) of [9]. The goal of the LCH 
is to capture features of quantum gravity in an effective 
theory which interpolates between General Relativity at 
low energies and a nonsingular spacetime at high curva- 
tures. The LCH has been studied in a variety of black 
hole and cosmological spacetimes [7] - [14] (also see [15] 
- [17]). The assumption of the existence of a limiting 
curvature is well motivated given the existence of a fun- 
damental length scale it (such as the Planck length i p i 
or the string length £ s ). If there is a fundamental length 
then it follows (from dimensional arguments) that all cur- 
vature invariants are bounded 



\R\ < i] 2 , 



\R aP R a0 \ < q 



\C aPlS C a ^ s \ < £ 
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(1.2) 
(1.3) 



In this letter we investigate the process of Hawking ra- 
diation of a nonsingular black hole in (1 + l)-dimensional 
dilaton gravity [11]. The construction involves concepts 
from the LCH. However, since the Ricci scalar is the only 
invariant in two dimensions, the role of the higher deriva- 
tive correction terms is facilitated by the dilaton and its 
potential. The resulting solution is nonsingular every- 
where and resembles the Schwarzschild black hole at large 
values of the radial coordinate. At r = the solution ap- 
proaches the de Sitter metric. Cosmologically desirable 
consequences may arise if the universe is generated at the 
interior of a black hole (see, e.g. [18] and the references 
therein) . 

Although it is natural to question the physical validity 
of results in 2D dilaton gravity, there are strong moti- 
vations for using the theory as a toy model of quantum 
gravity and to study the process of black hole evapora- 
tion. In particular, a certain limit of superstring the- 
ory yields a generalized dilaton theory as the effective 
action (see e.g. [28]). This simplified model allows for 
completely analytic calculations, possibly capturing non- 
perturbative features of the physics, evaporation and fi- 
nal state of realistic four-dimensional black holes. * Fur- 
ther motivation is provided by [29], where it is argued 
that the thermodynamics of black hole in D > 4 is effec- 
tively described by the thermodynamics of black holes in 
two-dimensional dilaton gravity theories. Specifically, it 
is shown that the Bekenstein-Hawking (BH) entropies of 
single-charged dilaton black holes and dilaton p-branes 
with an arbitrary dilaton coupling parameter in arbi- 
trary spacetime dimensions are exactly reproduced by 



*For some early references to two-dimensional gravity see 
[19] -[27]. 
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the BH entropy of the two-dimensional black hole in the 
associated 2D dilaton gravity model. We are also in- 
spired by the successful generalization of a nonsingular 
cosmological model in two dimensions [13] to a nonsin- 
gular four-dimensional model [7,8] using dilaton gravity 
and the limiting curvature construction mentioned above. 
Based on this encouraging evidence, we believe that the 
two-dimensional results presented in this paper should 
apply to more realistic four-dimensional models. 



II. A NONSINGULAR BLACK HOLE 

Here we reconstruct the nonsingular two-dimensional 
dilaton black hole solution presented in [11]. We 
begin with the most general Lagrangian for gravity 
with a scalar field (e.g., dilaton gravity) in (1 + 1)- 
dimensions [30]: 



C = V~9 (D(<p) R + Giv) (Vp) 2 + /%)) . 
By performing a conformal transformation 



and requiring that 



dp dip 



(2.4) 



(2.5) 



(2.6) 



it is possible to conformally transform away the kinetic 
term for ip and we may simply start with the action 

S = J d>x^—g (D(p) R + V(p)) , (2.7) 

where V(p) = e 2a ^H(p). t 

Varying this action with respect to ip and the metric 
tensor yields 



R 



dVjip) _ dD(ip) 
dip dp 
V(<p) g aP = 2(V 2 3q(3 - V a V p )D(p) , 



(2.8) 
(2.9) 



respectively. We may redefine tp so that D(tp) = 1/tp and 
assume a spherically symmetric and static metric 



ds 2 = -n(r)dt 2 +p(r)dr 2 . 



(2.10) 



In the "Schwarzschild gauge" p(r) = l/n(r) and the 
EOM become 



<p 3 V(<p) - Amp' 2 + n'p'p + 2nipip" = 



dV 

dp 

V(p)+p- 2 p'n' 



+ p- z n" = 
0. 



(2.11) 
(2.12) 
(2.13) 



Implementing the limiting curvature hypothesis is equiv- 
alent to identifying a class of potentials V(tp) so that R is 
bounded for all values of the radial coordinate from to 
oo. The form of the potential is restricted by considering 
the behavior of the system at low and high curvatures. 
We require n(r) to reduce to the Schwarzschild solution 
in regions of spacetime with small curvature. Mathemat- 
ically, 



n(r) — > (1 — 2m/r) 



when r — > oo, p — > . 

(2.14) 



From equation (2.8) it is easy to see that limiting V'(p) 
in the high curvature regime (i.e. r — > 0, ip — > oo) will 
result in the limitation of R. 

While there are many potentials that interpolate cor- 
rectly between the above asymptotic regimes, a simple 
potential which satisfies the above criteria is 



V(ip) 



2mA 3 ip 2 



1 + mA 3 l 2 ip 3 



(2.15) 



where A is a constant and I is a constant representing 
the limiting fundamental scale. * Using this potential in 
the EOM one finds the exact solutions 



(2.16) 
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where r$ is the location of the black hole horizon. The 
function n(r) is plotted below in Fig. (I). When the mass 
of the black hole is large, n reduces to (1 — 2m/ r) as 
expected. 



^Note that in [31], it is argued that under certain circum- 
stances conformally related theories may not be equivalent. 



'Note that this potential is only one possible potential that 
will force the curvature to be limited. While it is not a unique 
choice all potentials that satisfy the asymptotic conditions, 
limiting the curvature, will give the same result. For a more 
detailed explanation of the construction of V(p) we refer the 
reader to the original source [11]. 
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FIG. 1. A plot of n(r) for the values I = m = 1. (These 
values are used in the plots throughout the paper.) Far from 
the point r = the solution is indistinguishable from the 
Schwarzschild solution. The function n vanishes at the hori- 
zon r — ro = 2m. Unlike the Schwarzschild solution, n re- 
mains finite at r = 0. 

The Ricci scalar R (the only curvature invariant in 
two-dimensions) remains finite for all positive value of 
the radial coordinate r, 



R 



2m(2r 3 - l 2 m) 



(2.18) 



(r 3 + l 2 m) 2 

hence, the spacetime is singularity free (see Fig. (II)). 
Note that R changes sign at the value r — (l 2 m/2) 2 / 3 
and takes on its limiting value R m ax at r — 0, 
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FIG. 2. The curvature R remains finite for all r from r — 
to r = oo. 

Near r = the spacetime becomes de Sitter. § To see 
this, expand n around the point r = 0: 



§In the conventions of [11], the cosmological constant (A oc 
R) for de Sitter space is negative. 



n(r) 



-C+- + 0[(rf] 



(2.20) 



Ignoring terms of order r 3 and higher, the line element 
becomes 

ds 2 = -f^-c) dt 2 +(^-C\ dr 2 , (2.21) 

where C is a constant. Taking C = 1, we may define a 
new coordinate r given by 



dr 
dr 



(<7> 2 



1 



-1/2 



which gives 



I 



= arccosl 



(t) 



(2.22) 



(2.23) 



In the (t, r) coordinates the metric takes on the recog- 
nizable de Sitter form 



ds 2 = dr 2 - sink 2 - dt 



(y) 



(2.24) 



Another conceivably interesting region is near the fun- 
damental scale r ~ I. Expanding around r — I, n takes 
the form 



n(r) ~ C + 



2m 



1(1 + 1 



■(r-l)+0[{r-l) 2 }. (2.25) 



Here we choose a new coordinate x (and C = 1) defined 
by 



dx = , 2m(r-l) _ 1/2 
dr 1 1(1 + 1) ^ 



(2.26) 



1(1 + 1) 

In the (i, x) coordinates the metric becomes 

ds 2 = d X 2 - a 2 x 2 dt 2 , (2.27) 

where the constant a = m/l(l + 1). This is simply the 
Milne metric. 

One final comment is in order concerning the above so- 
lution. In [32], Myers and Horowitz argue that spacetime 
singularities should play a useful role in gravitational 
theories by eliminating unphysical solutions. Further- 
more, they suggest that any modified theory of gravity 
must have singularities in order to posses a stable ground 
state. The argument goes as follows: even if the theory 
claims to have a ground state with E < 0, one can always 
start with the Schwarzschild metric with m < E and ar- 
gue that it must be singular. Removing all singularities 
leads to the existence of states with arbitrarily negative 
energy. They suggest that a new mechanism must be 
found to prevent a state which resembles the negative 
mass Schwarzschild solution from existing in the theory. 

In the above model, we may argue that the negative 
mass Schwarzschild solution is not a well behaved solu- 
tion, in the sense that a singularity is present for m < 



3 



(see equation (2.18)). Furthermore, we may argue that 
the negative energy solutions may simply be thrown out 
due to their pathological, singular nature. Negative mass 
Schwarzschild solutions are not possible asymptotic solu- 
tions of the theory since a singularity is present at finite 
r. Thus, the instability of the theory for to < discussed 
in [32] does not occur in our theory. ** 



III. HAWKING RADIATION AND THE 
METHOD OF COMPLEX PATHS 

We now examine the process of Hawking radiation of 
the nonsingular black hole solution obtained in the previ- 
ous section. One of the simplest ways to do this is to em- 
ploy the method of complex path analysis introduced by 
Landau [33] and applied to Schwarzschild-like spacetimes 
by Srinivasan and Padmanabhan [34] (also see [35], for a 
related treatment). Here we use the method to calculate 
the Hawking temperature of a two-dimensional, static 
and spherically symmetric metric with horizon. We will 
then apply this formula to the black hole solution derived 
above. 

The line element is given by equation (2.10) in the 
Schwarzschild gauge p(r) = l/n(r): 



ds 2 = -n(r)dt 2 + n _1 (r)rfr 2 



(3.28) 



We assume that n(r) vanishes at some r and n'(r) is 
finite and nonzero at tq. This indicates that there is a 
horizon at the value r = tq. One must of course check 
that the singularity at r is only a mathematical one. In 
our case this is simple since the only curvature invariant 
in two-dimensions is R, which remains finite at rg (see 
equation (2.18)). Expanding n(r) around the point r 
gives 

n(r) = n'(r )(r - r ) + 0[(r - r ) 2 } = H(r )(r - r ) , 

(3.29) 

where we assume that H(ro) ^ 0. 

Now consider a scalar field which satisfies the Klein- 
Gordon equation 



□ 



mf 
h 2 



$ = 0. 



(3.30) 



If this field is propagating in the background spacetime 
given by equation (3.28), then 



1 <9 2 <f> d ( d$\ _ m,Q 

W + d~r { n{r) dr~ ) ~ ~W 



(3.31) 



The semiclassical wavefunctions satisfying the above are 
obtained by making the standard ansatz 



$(r,i) = e j fi 



(3.32) 



where S is a functional which will be expanded in powers 
of h. Substituting this ansatz into equation (3.31) gives, 



1 fdS\ 2 . . (dS^ 2 

-W)\to) +n(r) l^ 

_J_djS _ <PS_ _ dn(r) dS 
n(r) dt 2 Or 2 dr dr 



+ ml 



0. (3.33) 



We now expand S in a power series oih/i: 



S(r,t) = S (r,t)+*S 1 (r,t)+(jj S 2 (r,t) + . . . , 

(3.34) 

and substitute the result into equation (3.33). Neglecting 
terms of order j or higher we have 



n(r) 



dSo 
dt 



n{r) [w 



■m„ = . 



(3.35) 



This is simply the Hamilton-Jacobi equation satisfied by 
a particle of mass mo moving in a background spacetime 
with metric (3.28). The solution to this equation is 



So(r,t) = -£t± f Jt^yjst-rnlnir), (3. 



36) 



where £ is a constant and is identified with the energy. 
For simplicity we take too = 0. ft Now it is possible to 
solve equation (3.35) exactly. 

Using the usual saddle point method, the semiclassical 
propagator K(x",x') for a particle propagating from a 
spacetime point x" — (ti,r\) to a point x' = (t 2 , r 2 ) is 



K(x",x') = J\fe^ So{ - 



II I \ 
X ,X ) 



(3.37) 



where So is the action functional satisfying the classical 
Hamilton-Jacobi equation in the masslcss limit and J\f is 
a suitable normalization constant. With the solution to 
equation (3.35): 



S (x",x') = -£{t 2 -h)±£ f 



dr 
n{r) 



(3.38) 



we can calculate the amplitudes and probabilities of emis- 
sion and absorption through the event horizon at r$. 
Note that the sign ambiguity in equation (3.38) corre- 
sponds to the outgoing (<9<S > 0) or ingoing (dSo < 0) 
massless particles. 



**I would like to thank R. Myers for useful discussions con- 
cerning this point. 



ft Note that the results do not change for m ^ [34]. 
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If the points x" and x' are on opposite sides of the 
event horizon then the above integral diverges (since n _1 
diverges at r = ro). Therefore, to evaluate the integral we 
employ the calculus of residues and choose the contour 
over which the integral is to be performed around the 
point r = ro. 

Consider an outgoing particle at r — r\ < ro- The 
modulus squared of the amplitude for this particle to 
cross the horizon gives the probability of emission of the 
particle. Invoking the usual "ie" prescription, the con- 
tribution to iSo in the ranges (r, ro — e) and (ro + e, r^) 
is real. We take the contour to lie in the upper complex 
plane and find 

/"" Q+e dr 

So[emission\ = —£ hm / + real part 



/r -6 n ( r ) 



n(r ) 



real part , (3.39) 



where the minus sign in front of the integral corresponds 
to the initial condition that [dSo/dr) > at r = r% < ro, 
and lZ(ro) is given by equation (3.29). 

Now consider an ingoing particle with (dSo / dr) > at 
r = ri > ro- The modulus squared of the amplitude for 
this particle to cross the horizon gives the probability of 
absorption of the particle into the horizon. Choosing the 
contour to lie in the upper complex plane gives 



So [absorption] 



-£ lim 

e^O 



f r °~ € _dr_ 
' ro+e n(r) 



n(r ) 



- real part 



- real part . (3.40) 



This result agrees with the calculation of an outgoing 
particle (dSo/dr) > at r = r 2 > r . Here the contour 
is taken in the lower half-plane and the amplitude for the 
particle to cross the horizon is the same as that of ingoing 
particle due to time reversal invariance. 

Squaring the modulus to get the probability gives, 



P[ 



emission] cx e 



2ir£ 
hTZ 



(3.41) 



P[absorption] a e SK , (3.42) 

implying that 

47r£ 

P[emission] — e~ lnr P[absorption] . (3.43) 

Comparing this formula with the relation due to Hawking 
and Hartle [361: 



P[emission] = e^ P[absorptic 
we find the identification 

3 -i m\ 



4tt 



(3.44) 



(3.45) 



The method described above was shown to correctly re- 
produce the Hawking temperatures of the Schwarzschild 
black hole 



(3- 



87TTO 



(3.46) 



de Sitter spacetime and Rindler spacetime [34], and is 
valid for a variety of coordinate systems [37] . The method 
employs simple quantum mechanics and reproduces re- 
sults obtained by much more difficult means such as the 
calculation of Bogoliubov coefficients and other tradi- 
tional techniques, e.g. [36], [38], [39]. 



IV. AN ETERNALLY EVAPORATING BLACK 
HOLE 

Calculating the Hawking temperature of the black hole 
(given by the metric (3.28), with equation (2.17)) is now 
a simple task. Expanding n(r) near ro gives 



n(r) 



2mr . 

75 TV 

I m + rg 



r )+O[(r-r ) 2 



(4.47) 



In order to agree with the Schwarzschild solution we take 
ro = 2m, from which lZ(2m) — 4m/ (I 2 + 8m 2 ). Hence, 
the Hawking temperature of the black hole is given by 
(see equation (3.45)), 



T, 



km 



H 



ir(l 2 + 8m 2 ) ' 
This temperature is plotted in Fig. (3). 



(4.48) 



T H 




FIG. 3. The Hawking temperature Th asymptotically goes 
to zero as the mass goes to zero. 

The most important feature of this solution is that 
Th — * as m — > 0, in contrast to the standard 
Schwarzschild black hole which approaches infinite tem- 
perature in the zero mass limit (see equation (3.46)). No- 
tice that in the limit of large m (or in the limit as the 
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fundamental parameter I — > 0) the formula for Th appro- 
priately reduces to the Schwarzschild temperature. At 
the maximum temperature T max = h/A^/2itl, the curva- 
ture at the horizon, remains less than the magnitude of 
the limiting curvature. Hence, we are comforted that our 
semiclassical treatment of the radiation should remain a 
good approximation. 

A two dimensional version of Stefan's law gives the 
total power radiated by the black hole: 



dm 
~d7 



V T 



(4.49) 



The power is plotted in Fig. (IV). 



and VLHC). Most current discussions of black holes cre- 
ated in the lab are based on semi-classical calculations 
that are valid only when the mass of the black hole is 
much larger than M p i. When the black hole mass ap- 
proaches the fundamental scale the physics required to 
understand the process of Hawking radiation is rooted in 
a theory of quantum gravity. Our classical intuition con- 
cerning the creation of black holes in the lab may require 
refinement. (This is indicated by the above result). Fur- 
ther observational consequences of black hole remnants 
will be discussed in a future work [41]. 



MINIATURE BLACK HOLES AS DARK 
MATTER CANDIDATES 




FIG. 4. The power V ~ —dm/dt versus mass m. 



One of the greatest challenges of cosmology is to deter- 
mine the nature of the dark matter/energy which makes 
up most of the matter in the Universe. The eternal black 
holes we have described above could have very significant 
consequences for cosmology, since they have the poten- 
tial to contribute a sizeable amount of the dark matter 
in the Universe today (black hole remnants could behave 
as weakly interacting matter particles (WIMPS) [42]). §§ 
It is possible that in the early Universe geometric fluc- 
tuations produced a thermal (Boltzmann) distribution 
of black holes (see e.g., [46]). During the expansion of 
the Universe these black holes would slowly decay into 
radiation and very small black holes near or below the 
fundamental scale which would still be present today. 



Using the power it is possible to estimate the evapora- 
tion time of the black hole 



(P + 8m 2 ) 



2 i2 



\dm/dt\ 



(4.50) 



which is infinite in the limit m — > 0. In the limit of 
large mass (or taking I — > 0) the evaporation time re- 
duces to the familiar formula for the Schwarzschild black 
hole t, 



cvap 



m 3 . Our result implies that the black 
hole will radiate eternally. The mass will decrease as 
in the Schwarzschild case until most of the mass is radi- 
ated away, at which point the radiation decreases and 
one is left with a very slowly radiating, small black 
hole. H (At the maximum temperature the black hole 
has a Schwarzschild radius of r s = l/y/2.) 

Note that if this result applies to realistic black holes, 
then it should significantly affect the analysis of evaporat- 
ing black holes created at CERN's Large Hadron Collider 
(LHC) or in future collider experiments (such as CLIC 



^Note that a similar behavior was found in [40]. Here the 
authors considered a microcanonical treatment of black holes. 



VI. CONCLUSIONS 

In this paper we have explored the process of Hawk- 
ing radiation within the context of a specific model for a 
nonsingular black hole. We find that the black hole ra- 
diation mimics the Schwarzschild case for large values of 
the mass, but then reaches a maximum (determined by 
the fundamental scale) and then slowly decreases for all 
values of the time coordinate. The resulting miniature 
black holes could play an important role as dark matter 
candidates. 

A brief comment on the black hole information loss 
problem is warranted. The above model does not suf- 
fer from the traditional black hole information loss prob- 
lem. All the information which falls into the black hole 
is transported into the nonsingular black hole interior. 
In principle, an observer travelling into the black hole 
should be able to recover this information. Because the 



§§ For some early references on black holes as dark matter 
candidates see [43-45]. 
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black hole never completely evaporates, the information 
is never truly lost. 

A standard argument against the existence of stable 
black hole remnants is that it seems physically unrea- 
sonable for a large amount of information to be carried 
within a small (Planckian) volume [47] . One possible way 
to circumvent this issue was proposed within the context 
of charged dilatonic black holes [48] - [50]. Here it is 
suggested that what appears to be a small volume to an 
observer outside of the black hole is actually an infinite 
volume tube capable of storing infinite information. The 
scenario discussed in this paper shares this feature. The 
black hole is viewed from the outside as being a small 
remnant. The large interior core contains the missing in- 
formation which may be accessed only by travelling into 
the black hole. In [18] it was argued that the size of the 
universe inside the black hole is infinite. *** 

Of course, it remains to be shown that such physical 
properties will carry over to realistic four-dimensional 
black holes. This possibility is currently being investi- 
gated [51]. In order to generalize the model to four di- 
mensions we must use a higher derivative theory of grav- 
ity. The gravitational action should admit a solution 
which resembles the Schwarzschild black hole at large 
distances but with the singularity replaced by a de-Sitter 
universe. The method of construction will be similar to 
that in [8]. 

Our main conclusion is as follows. A specific construc- 
tion, based on the notion of a limiting curvature, is ca- 
pable of removing singularities in cosmological models. 
In this paper we have argued that this construction pre- 
dicts long-lived black hole relics with observational conse- 
quences and possibly solves the information loss problem. 

ACKNOWLEDGMENTS 

I would like to thank R. Brandcnbcrgcr, F. Lcblond, 
P. Martincau and R. Myers for enlightening discussions. 
I would also like to thank the Perimeter Institute in Wa- 
terloo, Canada for its hospitality during the completion 
of this work. 



***One problem which may remain in this scenario is that 
a large number of fundamental-mass particles will appear in 
loops and in the thermodynamic ensemble leading to a large 
degeneracy. It is conceivable, however, that they could be 
cured by suppressed amplitudes for creating such large cores 
if these amplitudes decrease sufficiently rapidly with the in- 
formation content [47]. 
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